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Abstract 

In this paper, we consider stochastic Schrodinger equations with two- 
dimensional white noise. Such equations are used to describe the evolution 
of an open quantum system undergoing a process of continuous measure- 
ment. Representations are obtained for solutions of such equations using a 
generalization to the stochastic case of the classical construction of Feyn- 
man path integrals over trajectories in the phase space. 
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Introduction 

In this paper we derive a representation for solutions of Schrodinger equa- 
tions with white-noise type coefficients (stochastic Schrodinger-Ito equations) 
using randomized Feynman path integrals over trajectories in the phase space 
(Hamiltonian Feynman integrals). These equations describe the limit behavior 
of a quantum system observed at discrete instants of time under the condition 
that the precision of measurement and the time intervals between measurements 
are proportional and tend to zero. Continuous observation of a quantum system 
can be defined (informally) as the limit of these observations. This enables one 
to assume that Schrodinger-Ito equations describe the evolution of a quantum 
system [20], [22] undergoing a process of continuous measurement. At the same 
time, they describe the so-called Markov approximation for the evolution of 
an open quantum systerrQ as opposed to the approximation given by Hudson 
Parthasarathy type quantum stochastic equations (see [12], [4], and [5]). 

An equation describing the evolution of a quantum system subject to a 
process of continuous measurement for some fixed observable (the operator of 
multiplication by a coordinate for an appropriate realization of the Hilbert state 
space in the form of i^(R^)) was first postulated in [25] to describe the sponta- 
neous reduction of the wave function. It was derived independently by Belavkin 

^By an open quantum system we mean one which is part of a larger quantum system. 
The evolution of such a system is described not by the Schrodinger equation but by a 
master-equation implied by it. Since master-equations are integro-differential, it is rather 
difficult to investigate them, and therefore use is made of various approximate equations. The 
Schrodinger-Ito equation can be regarded as one of these. 



[17] in the general situation. (In [25], use was made of the Hudson Parthasarathy 
quantum stochastic equations [26]; see also [18].) In the most important special 
case, this equation was derived independently by Diosi [21] at the physical level 
of rigor. For a derivation on the basis of the standard axioms of quantum me- 
chanics, see [8] (also [3], [15], and [16]). In the same paper [8] was announced 
the stochastic Schrbdinger equation with two-dimensional white noise. For a 
full derivation, see [32] . Apart from the local approach to the description of the 
behavior of a continuously observed quantum system resulting in an evolution 
equation that generalizes the Schrbdinger equation and takes into account the 
interaction between the quantum system and the measurement equipment and 
the effect of this equipment on the state of the system, there also exists a global 
approach developed in [30] and [31]. For the global description of the process 
of continuous measurement in the latter approach, a linear propagator for the 
quantum system is introduced in the form of an heuristic Feynman path integral 
over phase-space trajectories. 

For the relationship between these two approaches, sec [15] and [16]. Repre- 
sentations for solutions of stochastic Schrbdinger equations via Feynman path 
integrals were first obtained in [8], [2] and [16]. In these papers, the Feynman in- 
tegral was defined as the analytic continuation of the integral with respect to the 
Wiener measure (see [9]), as a result of which the analytic constraints imposed 
on the initial condition and on the potential turned out to be rather restrictive. 
Furthermore, a representation for the solution of the stochastic Schrbdinger 
equation with one-dimensional white noise was obtained in [13] under the con- 
dition that the potential and the initial condition in the Cauchy problem are 
the Fourier transforms of countably additive measures. In this approach, use is 
made of the definition of the Feynman path integral via Parseval's equation [7] , 
[14], [10], [9]. We use Feynman's original definition [11], [23] of the functional 
integral via the limit of finitely multiple integrals and extend the approach based 
on Chernoff's theorem [19] to the probabilistic case. This approach was first 
used in [35] to obtain a representation for the solution of the heat equation on 
a compact Riemannian manifold and then in [34] for the representation of the 
solution of the Schrodinger equation using a Feynman path integral over trajec- 
tories in the phase space. In the present paper, we obtain a representation for 
the solution of the stochastic Schrodinger equation using randomized Feynman 
path integrals over trajectories in the phase space (randomized Hamiltonian 
Feynman integrals). 

1 Pseudo-differential operators and the stochas- 
tic Schrodinger equation 

Definition 1. For an arbitrary r e[0, 1] we define a map A: H H from the 

space of complex-valued functions on M} x M} to the space of linear operators 
on 1,2 (K) as follows: the action of the operator H: D (h) c (R) (R) 
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on a function is given by the formula (see [34] and [24]) 
(^^) (<?) = ^ iTo/ / ^ " + e'^'°(«-«''V (90) dqodpo, (1) 

where the Umit is taken in (M). Let D (iij be the set of all ip € LP' (M) such 
that i-Lip exists. 

The function • , • ) is called the r -symbol (or the classical Hamiltonian 
on the phase space x M}) for the pseudo-differential operator %. The map 
A determines the r-quantization. 

The operator A for t= is called the operator of qp- quantization and, for t= 
1, the operator of pq- quantization. This terminology is related to the fact that 
the pseudo-differential operator corresponding to the T-symbol H{p, q) = pq is 
equal to qp for r= and to pq for t= 1. Here and henceforth, q and p are the 
coordinate and momentum operators given by the formulae 



q : f eDom{q)cL^iR)^[q^qf{q)], 
p : f€Dom(p)cL^{M)^[q^-if'{q)], 

respectively. 

The operation of r-quantization for r= 1/2 is called the Weyl quantization. 
We also note that if H{q,p) — f{q) -f g{p) for all p, 5 € M and some functions / 
and (7, then the result of quantization does not depend on the parameter r. 

In what follows, we consider the stochastic Schrodinger equation with two- 
dimensional white noise and interpret it as the Ito stochastic equation 



d^{t) = [(^-m-^^k^q)-^h\p)){^it)) 



dt 



- VT^Tfc (9) (t)) dWi (t) - yTJI^h ip) (^ (t)) dW2 (t) , (2) 

where H is the (internal) Hamiltonian obtained for the observed system by the 
r-quantization of the classical Hamiltonian H and where k{q) and /i(p)are the 
(non-commuting) differential operators corresponding to the rcal-valucd sym- 
bols {q,p) 7\k{q) and {q,p) h{p), respectively. Furthermore, Wi and 
W2 are independent standard Wiener processes and (t) G LP (M) is a random 
(wave) function describing the evolution of mixed states for the observed system. 
Equation (2) describes the evolution of an open quantum system undergoing a 
continuous measurement of the observables k{q) and h{p). It was considered in 
[8] and [32] in the special case when k{q) = g, h{p) = p and LI{q,p) = p^ /2+V{q) 
for all p,q & and some real-valued function V . 

For arbitrary functions h, k and LL, equation (2) can be obtained by the 
method used in [8] and [32]. For this, it sufBces to choose an appropriate 
realization of the Hilbert state space. 
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2 Feynman path integrals over trajectories in 
the phase space 

Let i? bo a real vector space, G <Z E* the space of linear functionals on E (it is 
assnmed that E and G satisfy the duality relation) and b a quadratic functional 
on G. Then, for all a G £" and a G C, the Feynman a-pseudo-measure with 
correlation functional h and mean value a is a generalized measure $6,0, « on E 
whose Fourier transform is given by the formula 



for all G G. 

Let Q and P be infinite-dimensional locally convex spaces which, as vector 

spaces, satisfy the relations Q = P* and P = Q* . li E = Q x P is the phase 
space and G = P x Q, then the zero-mean Feynman -pseudo- measure on E with 
correlation operator given by the formula b{p, q) = 2p{q) for all {p, q) G G is 
called the Hamiltonian Feynman integral (or the Feynman path integral over 
trajectories in the phase space). 

In this case, the value of the Feynman pseudo-measure on a given function 
/ : Q X P C (the Feynman integral of /) can be defined using the limit 
of finitely multiple integrals. Let {Qn}„ and {-Pn}„ be sequences of finite- 
dimensional subspaces (with dimQ„ =dimP„ = n) in Q and P, respectively. 
The sequential Feynman integral of the function / over trajectories in the phase 
space is defined as the limit (if it exists) of the expressions 



as n 00, where the integration is with respect to an arbitrary Lebesgue 
measure. 

For the representation of solutions of Schrodinger type equations by means of 
sequential Feynman integrals, special importance is attached to the case of the 
Feynman path integral over phase-space trajectories for some spaces Q and P of 
real- valued functions on the closed interval [0, t]. Let z G M'^, t g[0, 1] and t > 0. 
Let Q be the vector space of all real-valued functions on [0, t] whose generalized 
derivatives are measures on [0, t] and let P = {/ G C ([0, , M) : / (t) = 0}. In 
this case, the duality relation between Q and P is given by the formula 






for all G P and G Q, where (s) denotes the measure equal to the 
generalized derivative of the function (•). 
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Definition 2. The sequential Feynman integral 



JQxPO 



IQ 

of a function F : Q x P'^ C over trajectories in the phase space Q x is 
defined as the limit (if it exists) of the finitely multiple integrals 



In{F,z) = TT-^ / F{Jr{qo,--- ,qn),Jl{po,--- ,Pn)) 

{2-k) J^2^ 

X exp i ^ Pfe (9fe+i - <lk) \ dqodpo ■ ■ ■ dq„-idp„-i (3) 

as n — > oo. Hero p„ = 0, q„ = z and, for each r g[0, 1], the expression Jr 
is an (injectivc) map from M"+^ onto the space consisting of functions that 

are constant on each of the intervals (j^—fp^, fc = 1, • • • ,n. Furthermore, 

very n-tuple {qo,-..,qn) "= K""*"^, the function Jr{qo, qn) assumes the value 

(1 - T)qh + rqk-i on the interval (^^, f ) , A; = 1, • • • , n. 

Remark 1 (see [9] and [34]). For t= 0, the sequential Feynman integral 
defined above can be interpreted as an integral over the space of right-continuous 
functions. The case r= 1 corresponds to the Hamiltonian Feynman integral over 
the space of left-continuous functions and, for r g(0, 1), the set of functions that 
forms the convex hull (with weights (1 - t) and t) of the above function spaces 
should be regarded as the phase space. 



3 Feynman path integrals and representations 
of solutions of evolution equations 

The existence of Feynman path integrals over phase-space trajectories (that 
is, the convergence of the corresponding finitely multiple integrals) was proved 
in some special cases in [1] and [6] using the finite-difference method. The 
definitions of the Feynman integral via the analytic continuation of the integral 
with respect to the Wiener measure and via Parseval's equation were applied in 
[9] in the representation of solutions of Schrodinger equations. 

A new approach based on Chernoff 's theorem was first applied in [34] . This 
made it possible to extend substantially the area of application of Feynman's 
formulae (that is, the representation of solutions of Schrodinger equations using 
Feynman integrals). It was noted in [34] that if ip : IR+ >-)■ (M) is the solution 
of the Cauchy problem for the Schrodinger equation with initial datum ip^ and 
Hamiltonian equal to the pseudo-differential operator with r-symbol H, 



then the relation 



^(t)^e-**«^o- lini e-^^« (5) 

n— >-oo \ J 

is a representation for the solution of the Schrodinger equation via the Feynman 
path integral over trajectories in the phase space. Indeed, it can be verified that 
the right-hand side of (5) is a function whose value at a point z coincides with 
the limit of finitely multiple approximations of the Feynman integral 



1 exp { ^ {U («) ' («)) (0)) {^^1^ ■ 

The formula (5) was proved in [34] for a rather wide class of Hamiltonians using 
Chernoff 's theorem. 

In the present paper, we extend this approach to stochastic differential equa- 
tions of type (2). In this case, the corresponding Feynman formula changes. We 
shall show that if a random function ip : R+ i-^- (R) is the solution of the 
Cauchy problem for the stochastic equation 

dip = Aipdt + BipdW {t) (6) 

with initial datum p, where A and B are pseudo-differential operators on (R) 
and W is the standard Wiener process, then, under certain conditions on A and 
B, the relation 

ip{t) = lim hat fexp + BAWk,n + -a\\ipo (7) 

n-^oo \ I 2n n ) J 

holds, where AWfe,„ = W{tk/n) - W{t{k - l)/n) for all k, k = l,...,n, and 
hat(M) = M. The additional factors exp | — under the product sign cor- 
respond to Ito's formula. The right-hand side of (7) can be interpreted as a 
randomized Hamiltonian Feynman integral. 



4 A stochastic analogue of Feynman's formula 

The approach used in [34] to find solutions of non-stochastic Schrodinger equa- 
tions is based on the construction of a family of operators approximating in the 
sense of Chernoff (see [36]) the resolvent operator semigroup for the Schrodinger 
equation. Let Di be an essential domain for the operator "H, which means that 
the operator ^"H, D {iij ^ is the closure of {%, Dij . The one-parameter opera- 
tor family {S{t)}t>o approximates the semigroup with generator —iT-L (by defi- 
nition, it is precisely the resolvent semigroup of equation (2) for fi^ = jj^ = 0) 
in the sense of Chernoff if the relation S{t)f = f — itHf + o{t), t — ^ 0, holds for 
all f € Di. Then Chernoff 's theorem implies that the relation 
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e-''%^= lim (s(-\\ ^0 

n-^oo \ \n J J 

holds for alH > and ip^ G LP' (R). If e~'*'^ is taken as S{t), then we obtain 
the representation (5). 

This method cannot be used for the representation of solutions in the case of 
stochastic Schrodinger type equations (with /Ui,/i2 > in (2)) since the solution 
is a random function in (M) and the family of operators approximating the 
resolvent family for equation (2) is non-deterministic. We also note that, since 
the right-hand side involves a random process, the semigroup property does not 
hold for the resolvent family corresponding to the stochastic equation. Never- 
theless, if {T'r }s>r>o resolvent family of (random) operators on (M) 
that corresponds to the Cauchy problem for the stochastic equation (2) (that 
is, is defined by the formula T^ip^ = ip (s) for arbitrary s and s > r > 0, 
where (p is the solution of (2) with (p (r) = <^q, under the assumption that such 
a solution exists and is unique), then, since Wi and W2 are processes with inde- 
pendent increments, we have the following stochastic semigroup property: the 
distribution of depends only on the difference s — r. This enables us to gen- 
eralize the approach based on the notion of equivalence in the sense of Chernoff 
to the stochastic case. 

In what follows, we consider equation (2) for the case r= 0, that is, the 
operator H is obtained from "H by means of qp-quantization. It is assumed that 
the formula (5) holds for the Hamiltonian H. A sufficient condition [34] for the 
fulfillment of (5) is that the relation 

H{q,p) = ko{q) + ho{p)+l{q,p) (8) 

hold for all q,p € M and some real-valued functions ko,ho,l G (M). In 
addition, we assume that —il-L is the generator of a strongly continuous operator 
semigroup. 

We consider the family {Qr}o<r<s random operators on (M) defined 
by the formula 

(Q'rf) (9) = exp {-VMlfc (9) {Wi (s) - T^i (r)) - Mi (« - r) (9)} / (?) , g e M. 

(9) 

We claim that the function s ^ Q^f is the solution of the Cauchy problem for 
the equation 

it) = - yfc' (9) (*)) dt - VJhk {q) iv (*)) dWr {t) (10) 

with initial condition <^(r) = /. 

We take the total derivative of Q^f with respect to s, 

dsiQtfM = {-^/]lMqjdWi{s) + ^fi,k^q){dW^{s)f 
-Ml {s-r)k^ {q)ds}{Q:.f){q). 
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According to Ito's formula (dWi(s))^ = ds, it follows that s n- Q*/ is the 
solution of equation (10). The fulfillment of the initial condition is obvious. 

It can be similarly shown that if {-P^ }o<r<s operator family given by 

the formula 

(Prf) i<l) = ^ [ I exp {-VJ^2h iPo) {W2 (s) - W2 [r]) -fi2is- r) (po)} 
■^T^ Jr Jr 

then s H^- P// is the solution of the Cauchy problem for the equation 

(t) = -^h^ (P) m dt - (p) it)) dW2 (t) (11) 

with initial condition (f{r) = f. 

As shown in [34], if = e~'*'^ for all t, then the function t ^V'o approx- 
imates the solution of the equation 

dip {t) = -iHip {t) dt (12) 

in the sense of Chernoff. 

We note that the sum of the right-hand sides of the (linear) equations (10), 
(12) and (13) coincides with the right-hand side of (2), and it is therefore to 
be expected that the operator family {U^}^^^^^ determined by the formula 
Uf. = Q^Ys-rPr approximates, in a sense, the resolvent family for (2). 

The following relation will be called the Feynman stochastic formula: 

ToVo = - lim C/^_i),/„ • • • [/*/>o. (13) 

where w — lim denotes a kind of convergence (defined below) for (R)-valued 
random variables. To see an analogy with the usual Feynman formula, it suffices 
to note that t^T^Lp^ is the solution of equation (2) with initial datum ip^ and 
that the right-hand side of (14) is a finite-dimensional approximation to the 
Feynman integral of a random function (sec (7)). Finally, for fi^ = ij,2 = 0, we 
obtain the standard Feynman formula in [34]. 

Definition 3. Let be random variables defined on a probability space 
{Q,0,F) and having values in (R). A random variable ^ on the same space 
is the w-limit of ^„ as n ^ 00 if and only if IE jjC ~ CnllL2(R) tends to zero 
as n ^ 00. Here E denotes the operator of mathematical expectation in the 
probability space (fJ, G, P). 

To prove formula (14), it is necessary to consider asymptotic properties of 
the operator U^^ as .s — r — >• 0. In what follows, we use a well-known result of 
Levy [28] (see also [29]) on the local smoothness of a Wiener process. 
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Lemma 1. Let {W{r)}o<r<i be the standard Wiener process. Then the 
relation 

\W(r + v)-W(r)\ , 

hm sup — = 1 

"^"re[o,i] v—2vlnv 

holds with probability 1. 

It follows from Lemma 1 that the relation 

\W{r + v)-W{r)\ 
hm sup 1 = 

holds for each £ > with probability 1. 

2 

For g e M we write 62 (q) = e'^ — 1 — q— ^ . It is clear that there is a constant 

Ai> such that the inequality \eq {q)\< Ai\q\'^ holds as g — >■ 0. For brevity, in 
what follows we shall use the notation AWp"" = Wj {r + v) — Wj (r) for j = 1,2, 
and therefore the relation 



exp {-^fc (q) {Wi (r + v)- Wi (r)) - fi^vk^ {q)} 

= l-^k{q)AWr-f^ive (g) 

+ ^ [^.AWl'^k (q) + fx^ve (g)] ' + £2 {^^AWl^'k (g) + ii^v^ (g)) 

^ l-^,k (g) AWl^^-^i^vk^ (g)+i [^^AWl'^k (g) + ^i^vk^ (g)] \c {v, r, g) 

will hold for the function defining the action of the operator Q^^^ (see (9)). 

It is assumed below that the function k is bounded (sup^gg |A;(g) | = Kl < 00) 
and belongs to (M). We claim that the relation E \\c{v, r, •)IIi,2(m) = 
then holds uniformly with respect to r as ^ 0. By definition, 

c{v,r,q) = \nWk^{q) + {^^,f'\,AWl^^k^q) 
+£2 {^^AWl'^'k (g) + ii^vk^ (g)) . 
Therefore we have the inequalities 



Klv- 



1 



'2'KV 

2 



\z\e 



K?V-^y/v = 0(l 



3/2 



E ||£2 {^.AWl'^'k (g) + fi^vk^ (g)) | 



AiE 



£2 {^AWl^"k{q)+^i,ve (g)) 



< 



1,2 
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The expression on the right-hand side of the third of these inequaUties has 
an upper bound equal to the sum J2^=o^j^\^^i'^\'' 11^^^^"''^ IL2 

with some constants bj. Since EjAM^f'^p is proportional to v^^^ and we have 
||fcQ:||j^2 < 11^11^2 K^~^ for all a > 1, each of the terms in this sum is of order 
0(?;3/2). Consequently, E ||c (w, r, •)||^2 = 0{v^/^), that is, E \\c{v, r, ■)\\r(L2) = 
0(w3/2), where £ is the space of continuous linear operators in (M). The 
uniformity of the estimates in r follows from the fact that the distribution of 
AM^I^'^ does not depend on r. 

Everywhere below, with arbitrary s and r, s > r > 0, denotes the min- 
imal a-algebra relative to which the random variables — Wi(r)} are 
measurable. We note that, since Wiis a process with independent increments, 
the cr-algebras Qf.^ and Q^l are independent if the intervals (ri,si) and (r2,S2) 
are disjoint. 

Lemma 2. Let fc be a bounded function belonging to (K). Then for 
arbitrary r,v > 0, we have the following asymptotic expansion of the random 
operator QJ^^" with respect to the parameter v: 

Q;+^ = l~^,k {q) AWr - {q) + ^vk' (q) + O {v^"^) , v ^ 0, 

(14) 

where 0(w^/^) denotes a random operator in (M) measurable relative to Q^.'^^ 
such that the mathematical expectation of its norm is of order 0{v'^^^) uniformly 
with respect to r. Furthermore, (^^' " is a zero-mean random variable measurable 
relative to Q!^'^'" whose distribution does not depend on r or t;, and 1 denotes 
the identity operator in (M). 

Proof. It follows from the previous argument that 

Q;+- = 1 - VMlfc («) AT^r" - ii.vk^ + ^k' (AW-[''f + O (v^/^) 

= 1 - {q) AWl'^ - i,,vk^ + !^ve - ^ + o (v^/^) . 

It remains to note that = is a zero-mean random variable 

distributed according to the law (1) ~ 1- 

An asymptotic expansion for P^^" can be obtained in a similar way (ev- 
erywhere below, denotes the minimal cr-algebra relative to which the ran- 
dom variables {W2{t) — W2{r)}r<t<s are measurable). If h G (M) and 
supggg \ h{q)\ < 1, then the relation 

P;+^ = 1 - V^/i ip) AW^'^ - ^vh' (p) + ifvh' ip) rj^'^ + O {v'^') , v^O, 

(15) 

holds, where 0(w''/^) denotes a random operator in (R) measurable relative 
to V^~^'" such that the mathematical expectation of its norm is of order 0{v^^'^) 
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uniformly with respect to r. Furthermore, rj^'^ is a zero-mean random variable 
measurable relative to 7^^+" whose distribution does not depend on r or v. 

To prove (16), it suffices to note that the operator P^+'' can be written in 
the using the Fourier transform F in L^: 

P^+^ = Fexp {-^/lAVFa^'" - fX2vh'^} F~^- 

An expansion similar to (15) can be obtained for the exponential in this formula, 
and the application of the Fourier transform operator on the right and left of 
the exponential leads to the asymptotic expansion (16) since this operator in 
unitary and non-random. 

The asymptotic expansions (15) and (16), the definition of the Hamilto- 
nian symbol (8) and Lemma 1 imply an asymptotic expansion for the operator 

Jjr+v _ Qr+v^-ivUpr+v. 

[/;+" = 1 - VmI^ {q) AWT'" - Y^fc' (9) + Y«^' (?) C'" - VJ^^h (p) ^w^^" 

-^vh^ ip) + ^vh^ iP) r,'-'" - tv [ko (q) + ho (p) + /) 
+VMl7^fc (q) h ip) AWl'^AW'^'^ + o (v'/^-') , 

where e> is an arbitrary number. 

To derive a similar asymptotic expansion for the operator T^+^, we need the 
auxiliary assertion below. 

Lemma 3. Let A and B be operators in a Banach space X, let W denote the 
standard Wiener process and let (p : 1R+ i-^- X be a random function satisfying 
the Ito stochastic equation dip = Aipdt + B(pdW (t). Then for any £ > 0, the 
asymptotic expansion 



^{t) = ip{0)+B^{0)W{t)+Aip{0)t+^B^<p (0) (w {tf - +0 (ti-^) 
holds. Here o (^^~^^ denotes a random variable such that the relation 



t 0, 



lim = 

holds with probability 1. 

Proof. By the definition of a solution of the Ito stochastic differential equa- 
tion, the function (p satisfies the integral relation 

ip{t)-ip{0)= [ Aip{r)dr+ [ Bip{r)dW{r), (16) 
Jo Jo 

where /g Btp (r) dW{r) is the Ito stochastic integral. We shall find an asymptotic 
expansion for ip (t) as t ^ using the method of indeterminate coefficients. Let 

fit) = >p{0) + (3W{t) + at + ^7 (w {tf - + o (ti-'^ . 
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Wc substitute the right-hand side of this equation into the right- and left-hand 
sides of equation (18). Since we are interested in an asymptotic expansion to 

within o (t^~'^^, we can take into account only the constant in the asymptotic 

expansion obtained when ip (r) is substituted in J* Aip (r) dr (by Lemma 1, the 

other terms are of order o ). Consequently, 

J A(p{r)dr = A(p{0)t + o(ti~^Y 
We similarly conclude that 

B<p{r)dW{r) = Bip {0)W {t) + BpW {r) dW{r) + o (ti'"'^ 



Equating the left- and right-hand sides of (18), we arrive at the formula 



PW (t) + at + I [VF {tf - 1 
Aip (0) t + Bif (0) W (t) + is/? [W {tf 



1 



+ 



whence we obtain the relations a = A(p (0), /3 = B(p (0) and 7 = and the 
assertion of the lemma follows. 

Using Lemma 3 in the case of equation (2), we derive an asymptotic expan- 
sion for the operator T^+'' : 



nr+V 



+ ^vh^ ip) 1^^'^ - iv (fco (g) + (P) + i)+o (v^^^-^) 



{AW[-'')^-v 



where, as usual, = — ' and ri'^''" = — . We note that the di- 

rect application of Lemma 3 proves the expansion (19) only for r = 0. However, 
since Wi and W2 are Wiener processes, the distribution of the operator-valued 
random variable T^+'" does not depend on r, and the resulting expansion holds 
for an any r. Moreover, the expression o (t;^/^~^) on the right-hand side of (19) 
is an operator in (M) such that the mathematical expectation of its norm is 
of order o (t;^/^~^) uniformly with respect to r. 

As can be seen from the asymptotic expansions (19) and (17), the difference 
between T;+'" and is equal to the term y/JIuHk (?) h (p) AWl''"AW^''" + 

o (w'^/^^'^) which plays a key role in the proof of the randomized Fcynman for- 
mula (14). We also note that, in its standard form, the Chernoff theorem used 
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in the derivation of Fcynman type formulae involves, apart from the approxima- 
tion requirement for the operator family, the condition that the norm ||5'(i)|| of 
the operator in the approximating family should not exceed e'-'*, where C > 
is some constant (common to all values of t > 0). 

It turns out that a similar technical constraint is also necessary for proving 
the randomized Feynman formula. 

Lemma 4. Let h and k be bounded functions belonging to (M). Let the 
map t E IITqII be differentiable at zero. Then there is a constant C > such 
that, for all r, > 0, the estimates 

E||f/;+^|| < e<^'',E||T;+''|| < e*^^ 

hold, where||-|| denotes the norm on the space of linear operators in (R). 

Proof. It suffices to consider the case r = (the estimates for arbitrary 
values of r > will be the same as those forr = since the distributions 
of the operator-valued random variables under consideration depend only on 
the difference between the superscript and the subscript). We first derive an 
estimate for U^. If Ki = sup,gR |A; (g) | < oo, then \\Q^\\ < vT^TAM/i M _ 

Therefore E ||Q;i|| < e^^i'^i". 

Similarly, if A'2 = supg^jj 1^ (?) I < 00, then E ||Po^|| < es^^''^". To see this, 
it suffices to note that, by the definition of the quantization operator, we have 

Po" = F-^ exp {^h (•) W2 {v) - tx^vh^ (•)} F, 

where exp { y^ft (•) W2 (v) — iJ,2vh^ (•)} is the operator of multiplication by the 
corresponding function in (R) and F is the Fourier transform operator in 

(R). Then the inequality E||Po"|| < e^^^Ms" follows from the fact that the 
Fourier transform operator is unitary along with the estimate for the mathemat- 
ical expectation of the norm of the operator exp {^/Ji^h (•) W2 {v) — /ijf (')}■ 

We note that, since the operator- valued random variable Qq (Pq ) depends 
only on the realization of the process Wi {W2) respectively and the processes 
Wi and W2 are independent, we have the relation 



EIIC/o"!! <E||QSI 



-ivl-L 



EIIPo^ll < e^i'^i^e^^''^" 



We shall now prove a similar estimate for the operator Tq . Since {T'*}g<j,<^ 
is the resolvent family of random operators that corresponds to equation (2), 
the relation = T^T^ holds for all r e {0,v). The random variables and 
Tq are independent, and therefore the inequality 



EllTo^HEIlT^IIEIlToni 

holds. Since the distribution of depends only on v ~ r, the function / 
determined by the formula f{v — r) = E ||To || for all r and v, 0< r < v, is well 
defined. 

This inequality can be rewritten in terms of /: 
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.f{v)<f{r)f{v-r). 

In a similar way, the inequality f{v) < {f{v/n))" can be established for all 
V > and n e N. It follows from the hypotheses of the lemma that /'(O) exists 
and //(O) <1. Consequently, the relation 

f{v) < lim (/(i;/n))" = e-^'W^ 

n— foo 

holds. The lemma is proved. 

Theorem 1 (the randomized Fcynman formula). Let h and k be bounded 
functions belonging to (M). Let the map t i-^= E ||To || be differentiable at 
zero. Then, for an arbitrary function cpg S (M) and t > 0, the relation 



lim Ul 



-I jj{n-l)l/n TT>-l"-,„ 

(n-l)t/n'^(n-2)t/n ' ' ' ^0 V'o 



L/n 



holds. 

Proof. Everywhere below, f > is arbitrary and fixed and v = t/n. We 
write the difference between the operators on the left- and right-hand sides of 
(14) in the form 



JJ. 



(n— l)u 



(n-l)t>^(n-2)v 



j=n {j-l)v 



ttI Tj-kv 
J--lfe=j-l(^(fc_l)v- 



For all t; > and j = 1, n (of course, v and n are related by the formula 
nv = t), we define L'^ (R)-valued random variables by putting 



^3 



n 



j + l rpkv 

fe=n-'(fe-l)u" 



^0-1)^ TTl 



Then the proof of the randomized Feynman formula (14) becomes equivalent to 
the verification of the fact that 



E 



1 " 



as n — )• DO. 



(17) 



L2(R) 



To prove (20), it suffices to show that the expressions E ||^2^j^j are uni- 
formly bounded with respect to j in some neighbourhood of the point v = 
and that the limit relation E (Cj)Cfe)i2(R) — >■ as u ^ holds uniformly with 
respect to all j ^ k. For r,v > 0, let 3^^+" be the minimal |-algebra containing 
Q!:+'' and ^1+''. 
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Furthermore, for an arbitrary random variable ^, we denote by Ej;+"^ the 
conditional mathematical expectation E [^|3^o ^ • Thus, Ej;+" is the aver- 

aging operator over the cr-algebra 3^^+". 

Consider the expression 

r I T^'" - U^" \ 

llf:v\\2 _ / „j+l^kv 0-1)'' 0-1)" yrl rrkv l/N 

ll?jlll,2(R) - J^\'-h=n-^ik-l)v i^fc=j-l^(fe-l)v¥'0 I W 

X I n'''+^ T*^" 0-1)^ 0-1)" TTl rrkv („\ J„ 

It follows from the asymptotic expansions (19) and (17) that the relation 



^U-i)v ^0-1)" 1 



V 



holds. This relation and the estimates for the mathematical expectations of the 
operator norms in Lemma 4 imply the inequalities 



|2 



unUu) = ^K-i)Jk}i 

< e^''("-^)KM2 llfcll' WhfE^AWt'^^'^V^y 

< e^%,^,,\\k\\'\\hf + o{v'-'')]y,f. 

Since e> is arbitrary, the second of these inequalities provides the desired 
estimate. 

It remains to show that the covariances E (CJ) Cfe) j^2(]j) ^nd to zero as w — >• 
uniformly with respect to j ^ I. Consider the case j = n, Z = n — 1, for which 
the relation 




rp{n—l)v _ jj{n—l)v 

holds. 

We first find E"^^ j^^^ (Cn; Cn-i)^2(-j{-)i that is, we carry out averaging over 
the (T-algcbra 3^™_x)r- "^^'^ scalar product {^n^^n-i) contains only two random 
variables depending on J^(™Li)^, namely, ("-D"^ ("-D" ^nd T^^"_i^^- 
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li^p = ^l^j^iU^,^_-^^-j^, then the asymptotic expansions for these two random 
variables imply the relation 



- iv (fco (q) + ho ip) + i)+o (vi-^)] ^ ^"-^^'' V, 

IVJhIhk (q) h ip) Awt'^^'^Awt'^"'" + o (v'^- 



-{■n—l)v 



2)«V')l2(R)- 



To calculate E™_^^^ i^rn^n-i) L2(^gy it is necessary to remove the two pairs of 
square brackets in the above relation. It should be noted that the mathematical 
expectations of some of the terms resulting from the removal of the square 
brackets are zero. Therefore 



fvk^q)C 



■JT^h ip) AW^'^ - ^vh' (p) + ^vh' ip) T?'-'" 



-i{n—l)v 



u, 



{n—l)v 



-iv {ko {q) + ho ip) + i)+o {vi'^)] ^ ("-^)> ,o(.^-) c/;ri^»..(M) 



+ ((1 + 0(1)) 



p(n — l)v 
■ {n-2)v 



u, 



{n—l)v 
(n-2) 



rp(n-l)v _ jj{n—l)v 



+ 



It remains to note that < e*^*^" ||<^oll ^ 6*^* IIVolli whence it follows that 
the relation 



iE(c,e:-i) I = lEjEjrSEr^i). ia.c-i)} i 

<(l + o(l))o(«^-) e^^^ll^of 
+ {1 + (1)) o (vi-^) \\k\\ \\h\\ s/JHine^'^' W^of ^ 
holds as w — )• 0. 

The situation with arbitrary j and I is treated in a similar way. In this case, 
the convergence of E (C}",^;") to zero is uniform with respect to j ^ I since the 
asymptotic expansions for the operators T^j_i-^^ and J7q"_j)„ as w ^ do not 
depend on j. 
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5 Stochastic Feynman path integrals over tra- 
jectories in the phase space 

If a function on the phase space depends on a random parameter, then the 
Hamiltonian Feynman integral of that function can be defined in a natural way 
by extending the classical definition in [24]. Let {il,,Q,V) be a probability 
space. Everywhere below, the random Feynman path integral of a function 
F : C([0,t],R) X C([0,t],M) x O C such that F{£,p,xip,-) is measurable 
relative to (0,(7) for all S,g,S,p E C([0,t],K) over trajectories in the phase space 
is a random function in (K) equal to the w-limit as n — )• oo of the random 
functions z In {F, z) defined by equation (3). 

Theorem 2. Let h and k be bounded real-valued functions belonging to 
(R). Let ho and ko be functions from R to M, let / G (K) be a real- valued 
function and let 'H{q,p) = kO{q) + hO{p) + l{q,p) for all q,p £ M. Furthermore, 
it is assumed that if {r^}o<r<t resolvent operator family corresponding to 

equation (2) with Hamiltonian 'H obtained by the qp-quantization of H, then the 
map 1 E ||Tq|| is differentiable at zero. Under these conditions, the relation 

T^o = 1 exp |- ^* {tU (s) , ?p (s)) + (C, (s)) + ^i^h^ dp is))) dsj 

X exp |-V7^^* k is)) dWi is) - h is)) dW2 (.s)| 

x</'o(^,(0))$°'*'(rfe,,a!g (18) 

holds for an arbitrary function e (M) and all values of f > 0. 

Proof. By the definition of the stochastic Hamiltonian Feynman integral, the 

right-hand side of (21) is equal to w - lim„^oo Bl^n-i)/n^'t''(2-2'//2 ' ' ' ^o'^fo^ 
where for all r,v > the action of the operator on a function (f G (R) 

is defined as follows. Let 

fiq,p) = exp{- iiHiq,p) + fiik'^ iq) + ^12^^ ip))"" 

-^,k iq) AWl'^ -VJ^ fh ip) AW^'^} 

Jo 

for all q,p G K. Then B^+^c^ = fip, where A is the operation of qp-quantization. 

We claim that _B^+" coincides with Let /i iq,p) = exp{/i]^fc^ iq)v — 

y/jr[k iq) AW^-"}, /2 iq,p) = ex.p{-ivn iq,p)} and /a iq,p) = explii^h'^ ip) v - 
\/Jh lo ^ (p) ^^2 ^ } fo'^ q,p GR. Then, by the definitions of the correspond- 
ing operator families, we have the relations UJ^'^'" = /1/2/3 and = /i/2/3- 
The second of these was proved in [34] under the condition that /i, /2 and /a 
are non-random functions. The proof in the stochastic case is similar. 

For every functions' e i^(Mx (M)), we define an operator Jig) : (R) 1-^ 

(M) by putting 
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Then the qp-quantization of g can be written as 

g = J{g)F-\ 

The relations 

hhfz = J (/1/2/3) F-^ = A J ih) F-\h = /1/2/3 

follow from the fact that /i (/a) does not depend on p (g), respectively. 

Thus, the relation 5^;+" = C/^+'' holds for all r,v>0, and (21) follows from 
Theorem 1. 

Remark 2. The result in Theorem 2 can be extended to the case of r- 
quantization for an arbitrary t g[0, 1]. For this, it suffices to find asymptotic 
expansions as — >■ for the operator B^'^'" = / in the case when A is equal 
to the operator of r-quantization. It can be shown that these expansions will 
have the same form as (17). (Of course, in this case, the operators A on the 
right-hand side of this equation will correspond to r-quantization.) Then an 
analogue of Theorem 1 can be proved which will imply the representation of the 
sohition using the Hamiltonian Fcynman integral. The resulting formula for this 
representation will coincide with (21) except that the symbol is replaced 

by the pseudo-measure corresponding to the case of r-quantization. 

We also note that the use of methods based on Chcrnoff's theorem in the 
stochastic case also leads to randomized analogues of formulae in [9] , [33] . For 
this, suitable expansions into Dyson series must be considered for random mea- 
sures depending on realizations of the Wiener processes Wi and W2. 
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